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Variable-Mesh Method of Solving Differential Equations 


A multistep predictor-corrector method for the 
numerical solution of ordinary differential equations 
has been developed. The difference equations em- 
ployed are generalizations, for the case of variable 
mesh spacing, of previous formulas requiring fixed 
step size. In addition to retaining the high local ac- 
curacy and convergence properties of the earlier 
methods, the variable-mesh method was developed in 
a form conducive to the generation of effective criteria 
for the selection of subsequent step sizes in the step- 
by-step solution of differential equations. These cri- 
teria are based on considerations of truncation error, 
convergence of corrector iterations, and an extensive 
treatment of relative numerical stability. The algorithm 
is of importance in the solution of practical problems 
arising in engineering and quantitative scientific re- 
search. 

The variable-mesh multistep method has been tested 
by applying it to several single differential equations 
and to several systems of differential equations. This 
testing has given a fairly thorough demonstration of 
the effectiveness and reliability of the algorithm. One 
system of substantial importance for which the vari- 
able mesh approach proved especially effective was the 
problem of heat transfer to a supercritical fluid with 
variable physical properties and fully developed turbu- 
lent flow in a smooth tube. Another system was a sto- 
chastic model of enzymatically controlled cooperative 
unwinding and template replication of biological 
macromolecules. Most of the test .problems were se- 
lected because of their inherent potential, both in the 
behavior of the solutions and in the behavior of the 
partial derivatives of the right hand sides of the equa- 
tions with respect to the dependent variables, for pro- 
ducing numerical difficulties. Some are particularly 
suited to a variable-mesh treatment while others can 


be solved elliciently with constant-mesh increments. 
In the latter cases it is important to note that the ac- 
curacy obtained by the variable mesh method was 
about the same as that obtained using constant incre- 
ments with the same number of steps. This result indi- 
cates that the variable-mesh procedures do not have a 
degrading effect when they are used unnecessarily. 
Each equation was solved on the IBM System 360 
using single precision starting values and double pre- 
cision arithmetic to advance the solution. Values of 
5, the target ^ relative truncation error, ranging from 
10 to 10 were used for each equation. The ac- 
curacy obtained was roughly proportionate to the 
values of d specified. It was noted that the step lengths 
were limited almost entirely by the truncation error 
for the smaller values of 5 with the stabiliiy/conver- 
gence criterion becoming of increasing importance 
with increasing 5. 

Some of the problems were used in comparing the 
new algorithm with other fourth order numerical 
methods which also permit some variability in the 
mesh increments. The other methods used were the 
standard fourth order Runge-Kutta method, the 
Nordsieck method, and the basic Adams-Bashforth/ 
Ada ms- Moulton method, allowing doubling and 
halving of the increments with the latter. The new 
method proved superior to the other methods for 
these problems. 

Note: 

Documentation for the innovation is available 
from: 

Clearinghouse for Federal Scientific 
and Technical Information 
Springfield, Virginia 22151 
Price $3.00 

Reference: B69-10017 

(continued overleaf) 


This document was prepared under the sponsorship of the National 
Aeronautics and Space Administration. Neither the United States 
Government nor any person acting on behalf of the United States 


Government assumes any liability resulting from the use of the 
information contained in this document, or warrants that such use 
will be free from privately owned rights. 



Patent status : 

No patent action is contemplated by NASA. 

Source: R. VanWyk 
of North American Rockwell Corp. 

under contract to 
NASA Pasadena Office 
(N PO-10515) 
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